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Abstract 



Results are presented from a search for a narrow, spin-2 resonance decaying into a 
pair of Z bosons, with one Z-boson decaying into leptons (e + e~ or }i + u~) and the 
other into jets. An example of such a resonance is the Kaluza-Klein graviton, Gkk/ 
predicted in Randall-Sundrum models. The analysis is based on a 4.9 fb -1 sample of 
proton-proton collisions at a center-of-mass energy of 7 TeV, collected with the CMS 
detector at the LHC Kinematic and topological properties, including decay angular 
distributions as a novel feature of the analysis, are used to discriminate between sig- 
nal and background. No evidence for a resonance is observed, and upper limits on 
the production cross sections times branching fractions are set. In two models that 
predict Z-boson spin correlations in graviton decays, graviton masses are excluded 
lower than a value which varies between 610 and 945 GeV, depending on the model 
and the strength of the graviton couplings. 
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1 Introduction 

The standard model (SM) of particle physics does not include an explanation of the large dif- 
ference between the typical scales of gravity and the other fundamental forces (i.e., the "hier- 
archy problem"). This problem can be overcome by adding one warped spatial dimension to 
space-time, according to the Randall-Sundrum model (H El; in such a scenario, the existence 
of Kaluza-Klein excitations of a spin-2 boson, the graviton, is predicted. These models have 
two fundamental parameters: the mass Mi of the first graviton mode, referred to as Gkk in the 
following, and the dimensionless coupling strength k = k/ Mpi, where k is the curvature of the 
warped space and Mpi the reduced Planck mass (Mpi = Mpi/ \/8rz). 

In the Randall-Sundrum model discussed above (RSI), the graviton decay channels to dipho- 
tons and dileptons are the most sensitive, and several searches for such signatures of the RSI 
scenario have been performed at the Tevatron (3]-^J and at the Large Hadron Collider (LHC) 
(6]-|9|. However, different scenarios with warped extra dimensions (ED) allow the SM fields 
to propagate in the ED |Tl0l - [T2)| , resulting in a suppression of decays to diphotons and dilep- 
tons. These models are more compatible with electroweak precision tests and limits on flavor- 
changing neutral current processes than the original RSI. Furthermore, they are characterized 
by different couplings of the graviton to the SM fields. Such couplings result in two distinc- 
tive effects: the branching fraction to SM vector-boson pairs can become dominant for certain 
values of the model parameters, and an enhancement in the longitudinal polarization of the 
resulting vector boson pairs is predicted. 

The CDF and ATLAS experiments have searched for beyond-SM particles decaying to ZZ pairs 
in the 2£2q, 4£ and 2£2v final states fl~3l IT4| . without specific requirements on the spin or the 
production mechanism of the particle. This Letter presents a search for heavy narrow spin-2 
resonances in the ZZ final state and investigates both the original RSI model and the Agashe- 
Davoudiasl-Perez-Soni (ADPS) model [fTUI . The distinct angular distributions resulting from 
different polarizations in the decay are exploited to enhance the signal sensitivity. The process 
studied, Gkk — > ZZ — > qq£ + £~ (£ = e, combines a high branching fraction with good mass 
resolution and limited background rates. The dominant background is Z-boson production 
with associated jets ("Z+jets"). The cross section for this process falls rapidly with the invariant 
mass of the candidate Z-boson pair, mzz- 

2 Event reconstruction and selection 

We search for a fully reconstructed decay chain of the graviton Gkk — > ZZ — > qq£ + £~, where 
the charged leptons £ ± are either both muons or both electrons, and each quark is associated 
with a jet in the Compact Muon Solenoid (CMS) detector. This search covers a graviton mass 
range between 400 and 1200 GeV. At higher masses, the two jets from the hadronic Z-boson 
decay merge into one, therefore the efficiency of the present search signature is degraded. 

A detailed description of the CMS detector can be found in Ref. ||T5| . In the cylindrical co- 
ordinate system of CMS, (p is the azimuthal angle and the pseudorapidity (rj) is defined as 
rj = — ln[tan(0/2)], where 6 is the polar angle with respect to the counterclockwise-beam direc- 
tion. The central feature of the CMS detector is a 3.8 T superconducting solenoid of 6 m internal 
diameter. Within the field volume are the silicon tracker, the crystal electromagnetic calorimeter 
(ECAL), and the brass-scintillator hadron calorimeter (HCAL). The muon detectors are located 
outside the solenoid and are installed between the layers of the steel yoke of the flux return. In 
addition, CMS has extensive forward calorimetry in particular two steel /quartz-fiber hadron 
forward (HF) calorimeters, which cover the pseudorapidity range 2.9 < \rj\ < 5.2. 
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2 Event reconstruction and selection 



Although the main sources of background are estimated from data, Monte Carlo (MC) simula- 
tions are used to develop and validate the methods used in the analysis. Background samples 
are generated using either MadGraph 5 HH, SHERPA 1.13 CZl, or pythia 6.4.24 (18|. Signal 
events are generated at Leading Order (LO), for both the RSI and the ADPS model, using a 
dedicated generator described in Ref. ||T9ll . Parton distribution functions (PDFs) are modeled 
using the parameterizations CTlO at next-to-LO |20| and CTEQ6 [l2"T| at LO. For both signal and 
background MC, events are simulated using a GEANT4-based model [22] of the CMS detector 
and processed using the same reconstruction algorithms as for data. 

Muons are measured with the silicon tracker and the muon system [23|. Electrons are identified 
as tracks in the tracker pointing to energy clusters in the ECAL |24| . Muons and electrons 
are required to have the component of momentum transverse to the pp beam direction, pj, 
greater than 20 GeV. At least one lepton must have pj > 40 GeV. Leptons are measured in 
the pseudorapidity range \t]\ < 2.4 for muons, and \t]\ < 2.5 for electrons. Electrons in the 
transition range between the barrel and endcap, 1.44 < \tj\ < 1.57, are excluded. Both the 
pj and r\ requirements are consistent with those in the online trigger selection requiring two 
charged leptons, either electrons or muons. 

Muons are required to be isolated from hadronic activity in the detector by restricting the sum 
of transverse momentum (energy) in the tracker (ECAL and HCAL), within a surrounding cone 
of AR = v / (Arj) 2 + (Acp) 2 < 0.3 (excluding the muon candidate itself), to be less than 15% of 
the measured pj of the muon, where At] and Acp are the differences in pseudorapidity and in 
azimuthal angle. Electron isolation requirements are similar, but vary depending on the shape 
of the electron shower. The details of electron and muon identification criteria are described in 
Ref. H2H1- 

Jets are reconstructed with the particle flow (PF) [26j algorithm. Reconstructed particle can- 
didates are clustered to form jets with the anti-fcx algorithm [27] with a distance parameter of 
0.5. Jets are required to be inside the tracker acceptance (\t]\ < 2.4) and to have pj > 30 GeV. 
Jets that overlap with isolated leptons within AR = 0.5 are removed. Jet energy corrections 
are applied to account for the non-linearity due to instrumental effects. These corrections are 
based on in-situ measurements using dijet, 7+jet and Z+jets data samples [28]. A quality selec- 
tion, primarily based on the energy balance between charged and neutral hadrons, is applied to 
remove spurious jets due to detector artifacts. The jet four-momentum is corrected for the con- 
tribution of overlapping minimum bias events coming from different proton-proton collisions 
(pileup) based on the effective jet area, following Ref. p9| . 

Each pair of oppositely charged leptons and each pair of jets passing the above selections are 
considered as Z candidates. Background suppression is primarily based on the dilepton and 
dijet invariant masses, mn and m». The requirement 75 < rrijj < 105 GeV is applied in order 
to reduce the Z+jets background, and 70 < tnu < 110 GeV to reduce background without a 
genuine Z decaying to leptons, such as tt. 

The kinematics of the jets is corrected by constraining the invariant mass of the dijet system 
to the mass of the Z boson. Experimental resolutions, typically lower than 10%, are taken 
into account in the kinematic fit. The statistical analysis is based on the invariant mass of the 
graviton candidate, rtizz, calculated by imposing the kinematic constraint to the dijet system. 

2.1 Jet flavor analysis 

As the dominant background is Z+jets (Z — > l + £~) production, we use the flavor content of 
the jets compared to those from Z-boson decays to discriminate against it. In Z+jets events the 



2.2 Angular analysis 
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jets will predominantly be produced by the fragmentation of low-mass quarks and gluons. As 
a result, we can exploit the fact that, in ZZ events with one Z-boson decaying into jets, there 
will be a larger fraction of jets from heavy-quark fragmentation than from Z+jets, and there 
will be no jets from gluon hadronization. We take advantage of both features by tagging the b 
flavor and introducing a likelihood discriminant that separates gluon and light-quark jets on a 
statistical basis, as described in Ref. II30II . 

To identify jets originating from the hadronization of bottom quarks, we use the CMS track 
counting high-efficiency (TCHE) b-tagging algorithm |3T] |32|| . The data are split into three 
b-tag categories. The 2 b-tag category is required to have one jet with a value of the TCHE 
discriminant passing the medium working point selection (~65% efficiency) and the other jet 
passing the loose working point selection (~80% efficiency). Events not selected in the 2 b- 
tag region are categorized as 1 b-tag if they have at least one jet satisfying the loose b-tagging 
requirement. The b-tag category contains all the remaining events. 

For the b-tag category, which is dominated by the Z+jets backgrounds, a quark-gluon discrim- 
inant (QG) is constructed from observables related to the multiplicity and distribution of the 
jet constituents |30|. The QG-discriminant values range from (gluon-like) to 1 (quark-like). 
We apply a loose requirement, QG > 0.1, in the b-tag category only. 

In order to suppress the substantial tt background in the 2 b-tag category, a discriminant, A, 
is constructed using resolution functions of the PF jets as described in Ref. II33II , to provide 
a measure of the missing transverse energy (#x) significance. We apply a loose requirement, 
2 In A(#t) < 10, in the 2 b-tag category only. 

2.2 Angular analysis 

Due to the tensor nature of the graviton, the angular distribution of its decay products shows 
a distinct signature. Five angles [6*, Q\, 62, Oi, O) defined in Ref. [19| fully describe the kine- 
matics of the qq, gg — > Gkk — > ZZ — > qq£ + £~ process. Further kinematic selection exploits 
these five angular observables, which are only weakly correlated with the invariant masses of 
Gkk and the two Z bosons, and with the longitudinal and transverse momenta of the graviton 
candidate. 

We construct an angular likelihood discriminant (LD) based on the probability ratio of the 
signal and background hypotheses, as described in Ref. [30J. These probabilities are param- 
eterized as functions of nizz- the full five-dimensional theoretical distributions multiplied by 
acceptance functions are used for the two signal types, while empirical functions describe the 
background probabilities. The left panel of Figure [I] shows the distribution of cos 9* after the 
pre-selection requirements, where 9* is the polar angle of one of the Z boson momentum vec- 
tors in the ZZ rest frame, measured with respect to the beam axis. Because the distributions of 
the signal and background events are so different, cos 9* is an important part of LD. The right 
panel of Figure [l] presents the LD distribution after the same pre-selection requirements.The 
expectation of longitudinal polarization in the ADPS model results in a visible enhancement of 
the discriminating power of the LD with respect to the Z+jets background. 

The figure of merit proposed in Ref. [34] is used to optimize the selection on the angular like- 
lihood discriminant in order to maximize the expected sensitivity to the production of an RSI 
graviton. Because of the different background contributions, the optimization is performed 
independently for the three b-tag categories of events. The results are approximately indepen- 
dent of the graviton mass hypothesis and lead to cut values of LD > 0.28, 0.35, and 0.21 for the 
b-tag, 1 b-tag and 2 b-tag categories, respectively. The optimization is performed only on RSI 
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2 Event reconstruction and selection 




Figure 1: Distributions of the cos 6* variable (left) and the angular likelihood discriminant, LD 
(right), after pre-selection requirements. Points with error bars show the data distributions. 
Stacked, filled histograms show the background expectations from simulated events with the 
different component breakdown. The contribution from Z+jets dominates the background. 
Solid- or dashed-line histograms (overlaid) show the expected distributions for RSI and ADPS 
gravitons with mass M\ = 700 GeV and k = 0.10, scaled up by arbitrary factors for illustration. 

signal and the same selection criteria are used for both models under scrutiny. 

Table 1: Summary of optimized kinematic and topological selection requirements. 

Preselection 

Pt(^) p T > 40 (20) GeV 
px(jets) > 30 GeV 
M(0 e ± : \t]\ < 2.5, ^ : \rj\ < 2.4 
| //| (jets) < 2.4 



b-tag 


1 b-tag 


2 b-tags 


b-tag 


None 


One loose 


Medium & loose 


Ang. LD 


> 0.28 


> 0.35 


> 0.21 


QG discr. 


> 0.10 


none 


none 


21nA(# T ) 


none 


none 


< 10 


nijj 




[75, 105] 


GeV 



m u [70, 110] GeV 

m zz [375, 1500] GeV 



2.3 Summary of selection 

The selection requirements are summarized in Table [l] When an event contains multiple Gkk 
candidates passing these selection requirements, the one with the largest number of b-tagged 
jets is retained, thus giving preference to the purest category. Further ambiguity between mul- 
tiple candidates is resolved by selecting the candidate with mu (before the kinematic fit) and 
mu values closest to the Z boson mass mz- Signal selection efficiencies depend on the graviton 
mass considered and range approximately from 3 to 10% (3 to 21%) for the RSI (ADPS) model, 
respectively. 
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3 Event analysis 

3.1 Background estimation 

In order to minimize systematic uncertainties associated with the background shapes, we es- 
timate the background distribution from the data. The mzz distribution of the selected events 
is split into three samples based on the number of b-tagged jets. The selected events are then 
examined comparing to 33 equally-spaced hypothetical graviton masses in the range between 
400 GeV and 1200 GeV. 

A signal window in the mzz distribution is defined around the expected position of the signal 
peak. The signal shape is parameterized using the sum of a double Crystal-Ball function (i.e., 
a Gaussian distribution with power-law tails on both sides) [35J and an empirical line-shape 
derived from the observed distribution of events with misassigned jets (a triangular shape con- 
voluted with a Crystal-Ball function). The latter component describes both mismeasurements 
and incorrect reconstruction of the graviton decay products. The size of the signal window 
corresponds to three times the standard deviation of the core distribution of the signal shape, 
which ranges from 13 (at Mi = 400 GeV) to 38 GeV (at Mi = 1200 GeV). For every mass value 
considered, the mzz distribution is fitted in the [375, 1500] GeV range after excluding the signal 
window, in order to estimate the background expectation from the sideband regions. The dis- 
tributions of mzz for the background are parameterized with exponential functions. To gain 
stability in the fits, we combine the electron and muon samples, assuming that the background 
shapes are the same. 

Figure |2]shows the result of the fit for the specific case of mzz = 700 GeV, overlaid with the mzz 
distributions split by both b-tag category and leptonic channel. While the normalization of the 
fit is determined independently for each plot, the good agreement between data and fit justifies 
the use of the same fit function in the eejj and the fijijj channels. The MC prediction normalized 
to the integrated luminosity of the data sample is also included in the plot. Although the agree- 
ment of the MC prediction with the data is fair, differences are observed in both normalization 
and shape. The most significant discrepancy between data and MC, in the b-tag category 
(11% in the total yield), is because the relative fraction of quark- and gluon-originated jets is 
not well reproduced in MadGraph MC, which affects the QG selection efficiency. Because the 
background is determined directly from the sidebands in the data, the results are not sensitive 
to the limitations of the simulated samples. 

The yields predicted by the background estimation from control samples in data are compared 
to the observed ones in Table [2] Additionally, the expected number of events for the hypothet- 
ical graviton signals in the RSI and ADPS models are listed for some benchmark mass values. 

3.2 Systematic uncertainties 

The uncertainties on the normalization of the background fit are statistical in nature since they 
depend on how many events populate the sidebands. The shape uncertainty is determined 
by propagating the uncertainty on the slope parameter of the exponential into the background 
estimation. Because of the steep fall of the mzz spectrum, these uncertainties are larger for low 
values of the graviton masses. In the case of the search for a 400 GeV (700 GeV) resonance, the 
uncertainties on the normalization and shape of the background are 10% (4%) for the b-tag, 
10% (5%) for the 1 b-tag, and 30% (15%) for the 2 b-tag category. 

The systematic uncertainties on the signal normalization are summarized in Table [3] We con- 
sider effects from lepton energy scale, resolution, selection, and trigger; jet efficiency, energy 
scale and resolution after calibration; pileup; significance requirements; heavy-quark flavor 
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3 Event analysis 
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Figure 2: The mzz invariant mass distribution after final selection in three categories: b- 
tag (top), 1 b-tag (middle), 2 b-tag (bottom). The two channels eejj and f/f//; are presented 
separately on the left and on the right, respectively. Points with error bars show distributions 
of data, stacked filled histograms show the background expectation from simulated events 
with the different components illustrated. Also shown is a hypothetical RSI graviton signal 
with mass of 700 GeV and k = 0.10. The red (lighter) solid line shows the central value of 
the prediction of background from the sideband extrapolation procedure (excluding a signal 
window of ±57 GeV around 700 GeV in this example). 
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Table 2: Observed and expected yields in the 4.9 fb~ dataset at y/s = 7TeV. The yields are 
quoted in the range 375 < mzz < 1500 GeV (excluding a signal window of ±57 GeV around 
700 GeV in this example). The signal expectations for the RSI and ADPS models are reported 
for three benchmark M\ points. The uncertainties shown are the sum in quadrature of the 
systematic and statistical uncertainties. 



b-tag 


1 b-tag 


2 b-tags 


Obs. yield 658 
Exp. bkgnd 652 ±25 


542 
548 ±23 


51 
51±7 


Signal expectation (RSI graviton, £=0.10) 


500 GeV 268±24 
700 GeV 50±5 
900 GeV 8.4±1.2 


212±16 
39 ±4 
7.2±0.8 


100±12 
12.4±1.6 
1.36±0.24 


Signal expectation (ADPS graviton, £=0.10) 


500 GeV 120±10 
700 GeV 11.3±1.2 
900 GeV 0.78±0.10 


102±6 
10.9±0.7 
0.71±0.05 


42±5 
3.2±0.4 
0.097±0.012 



tagging and quark-gluon discrimination; graviton production mechanism; and the luminosity 
measurement. 

Reconstruction efficiencies for leptons and their uncertainties are evaluated from control sam- 
ples in data [25]. The systematic uncertainties on jet reconstruction are evaluated by variation 
of the jet energy scale and resolution within calibration uncertainties. Uncertainties due to 
pileup are taken as the difference between reconstruction and selection efficiencies with the 
number of pileup vertices shifted by one unit up and down with respect to the average value 
measured in data. The requirement on the discriminant A(#j) translates into 1-2% reduction 
in signal and the resulting uncertainty is taken as the full inefficiency value. Uncertainties on 
the quark-gluon selection efficiency are evaluated using a selected 7+jet sample in data, which 
predominantly contains quark jets. The uncertainty on the b tagging has several sources, such 
as pileup effects and dependence on the fraction of b production from gluon splitting, and has 
been evaluated on a sample of dijet events enriched in heavy-flavor by requiring a muon to be 
associated with one jet II32I . 

For a given production model, there is an additional uncertainty in the production mechanism 
coming from the choice of PDFs. PDFs affect the distribution of the longitudinal momentum 
of the graviton, hence the signal acceptance. We follow the PDF4LHC 11201 136T139I recommen- 
dation to estimate the uncertainty due to PDF knowledge and to calculate the uncertainty on 
signal acceptance. The relative uncertainty on the integrated luminosity [40J is also applied to 
the signal normalization. 



4 Results 

Based on the normalization and shape of the mzz distribution for signal and background, ob- 
tained as described in Section 3.1 we perform an unbinned statistical analysis of the results 
using the same formalism as in |30ll4"T| . Since no significant excess is observed, we set exclu- 
sion limits at the 95% confidence level (CL) on the product of the graviton cross section and the 
branching fraction of Gkk ZZ. The upper limits are calculated using the CL S ||42ll43 l method 
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4 Results 



Table 3: Summary of sources of systematic uncertainties on signal normalization. All are ex- 
pressed as relative uncertainties on the cross section. 



Source 


b-tag 


1 b-tag 


2 b-tags 


}i reconst. & trig. 




2.7% 




e reconst. & trig. 




4.5% 




Jet reconst. 




2-4% 




Pileup 




1-4% 




#T significance 






1-2% 


b tagging 


2-7% 


3-5% 


10-11% 


QG-discrimination 


4.6% 






Acceptance (PDF) 




~1% 




Luminosity 




2.2% 





and compared to the theoretical cross sections times branching fraction of Gkk —> ZZ. 



Figure [3] (left) shows the upper limits for the RSI model and compares them to LO theoretical 
expectations for two benchmark values of the coupling strength k, equal to 0.05 and 0.10. These 
results allow us to exclude at the 95% CL gravitons with M\ < 945 GeV for k = 0.10, improving 
on previous analyses [14J where a phase-space decay model of the graviton is used. For a value 
of the coupling strength k = 0.05, we exclude the separate ranges Mi < 720 GeV and 760 < 
Mi < 850 GeV, because of an under-fluctuation of the observed limit around Mi = 800 GeV. 
The upper limit for a signal generated according to the ADPS model is shown in Fig. [3] (right). 
The limit on the production cross sections times branching fraction of Gkk — > ZZ, as a function 
of the graviton mass, is compared to the model predictions at LO, resulting in exclusion ranges 
of Mi < 720 GeV and Mi < 610 GeV for k = 0.10 and 0.05, respectively. 
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Figure 3: Observed (black solid lines) and expected (black dashed line) 95% CL upper limits on 
the cross section times branching fraction for Gkk ZZ, obtained with the CL S technique. The 
graviton signal was generated according to the model in [1J (RSI, left) and |fT0| (ADPS, right). 
The 68% and 95% ranges of expectation for the background-only model are also shown with 
green (darker) and yellow (lighter) bands surrounding the expected upper limit line, respec- 
tively. The predicted product of the graviton cross section and the branching fraction is shown 
as a red solid (dashed) curve for k = 0.10 (k = 0.05). 
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5 Summary 

Results have been presented of a search for a KK graviton decaying into ZZ, with one Z-boson 
decaying leptonically the other hadronically. CMS data at ■s/s = 7 TeV have been used, cor- 
responding to a total integrated luminosity of 4.9 fb . As a novel feature of this study the 
analysis targets specific graviton models, by exploiting information on angular distributions of 
the decay products to build likelihood-ratio discriminants that enhance the signal sensitivity. 

No excesses in the mzz distributions over the expected SM backgrounds are found, for the 
range 400-1200 GeV. Upper limits on the production cross section times branching fraction as 
a function of graviton mass are set at 95% confidence level. Exclusion ranges for the RSI model 
are Mi < 945 GeV for k = 0.10, and Mi < 720 GeV and 760 < Mi < 850 GeV for k = 0.05. 

For the first time, upper limits are obtained specifically for the ADPS model, which is more 
compatible with current indirect experimental limits than the original RS model. We place 
exclusion ranges of Mi < 720 GeV and Mi < 610 GeV for k = 0.10 and 0.05, respectively. 
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